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Abstract 

Concerning state estimation, we will compare two cases. In one case we cannot use the 
quantum correlations between samples. In the other case, we can use them. In addition, 
under the later case, we will propose a method which simultaneously measures the complex 
amplitude and the expected photon number for the displaced thermal states. 

1 Introduction 

Quantum estimation is essentially different from classical estimation regarding the following 
two points. The first point is that we cannot simultaneously construct the optimal estimators 
corresponding to respective parameters because of non-commutativity between them. It has 
been a serious problem since the beginning of the quantum estimation [Q, ^, |3). The second 
point is that we can reduce the estimation error under the assumption that we can prepare 
independent and identical samples of the unknown quantum state. It was pointed by Nagaoka 
Q H concerning the large deviation theory in one-parameter estimation. The purpose of this 
paper is to clear the second point concerning the mean square error (MSE). 

Our situation is divided into the following two cases. In the first case, we estimate the 
unknown state by independently measuring every sample. In this case, we may decide the 
n-th POVM from n — 1 data which have been already given. In the second case, we estimate 
the unknown state by regarding n-sample system as a single composite system. In this case, 
we may use POVMs which are indivisible into every sample system. In order to construct 
these POVMs, we need to use quantum correlations between every sample. The former is 
called the non-quantum correlation case and the later the quantum correlation case. When 
the unknown state is a pure state, the errors of both are asymptotically equivalent in the 
first order Q. Concerning the spin 1/2 system, see Hayashi J7j. 

In this paper, we formulate a general theory for the asymptotic quantum estimation. It 
is applied to the simultaneous estimation of the expected photon number and the complex 
amplitude for the quantum displaced thermal state. 

2 Asymptotic Estimation Theory 

In this paper, we use a quantum state family S parameterized by finite parameters 9 1 , • • • , 6 d : 

S := { Pe G S{H)\9 = (9\...,9 d )eec H d }, 

where the set S(H) denotes the set of densities on TL. For simplicity, we assume that pq is 
nondegenerate. 
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2.1 Non-quantum correlation case 

The non-quantum correlation case is formulated as follows. A pair £ n = ({-Mfe}? =1 , 9 n ) is 
called a recursive estimator where 9 n is a function estimating the unknown parameter from 
n data, and {Mk}^ =1 is a sequence of POVMs M%, il^u^i), . . . , M n (uJi, . . . , us n -\) as follows: 
the n-th POVM Mk(ux, . . . , Wfc-i) is determined by k — 1 data which have been already given. 
A sequence {£ n }^ =1 of recursive estimators is called a recursive MSE consistent estimator if 



lim 

n— >oo 



P£" ( dwi, «4j 2) • • • , diO n ) = 0, V6< 6 0, 



where 

P p n (dux, dio 2 , ... , duj n ) = tr pMx(du;i)tT pM 2 (uji)(du2) ■ ■ .tvpM n (ui, . . . , w„_i)( dw n ). 

We define the non- quantum- correlational Cramer-Rao type bound Cg (G) for a weighted 
matrix G(G is a <i x d real positive symmetric matrix.) as: 

C^® C (G) := inf jliminf Ti:nGVg(£ n ) {£ n }^ =1 is a recursive MSE consistent estimator j , 
where the MSE matrix V$(£ n ) is given by: 



We have the following equation: 



C, 



NQC 



(G) = mil TrG J ( 



7iU 



M is a POVM on W 



where J^ 1 denotes the Fisher information matrix at 9 of {tr pgM( dw)\9 € 0}. It is derived by 
Jensen's inequality H|. Under some regular condition, we show that there exists a recursive 
MSE consistent estimator {£ n }^ =1 such that ||: 



nTrGV e {£ n ) -» Cf QC (G) as n — > oo, V# G 9. 

According HolevoQ, we have Ji^ < Jg, where Jg is the Fisher information matrix 

defined as: Jg-i : j : = tr(Lg-i)* pgLg-j, Lg-i := (pg)^ 1 ^- Therefore we have the following 
inequality 



C^ QC {G) > Cf (G), 



where 



Cq"{G) := inf |trGV|V is a d x d real symmetric matrix V > J g 1 j 



TrGRe J7 1 +Tr 



G Im J. 



(9 
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2.2 Quantum correlation case 

Next, we formulate the quantum correlation case. For this purpose, we consider a quantum 
counterpart of independent and identically distributed condition. If Hi, ■■ ■ , H n are n Hilbert 
spaces which correspond to the physical systems, then their composite system is represented 
by the tensor Hilbert space. 

H W :=Hi®---®H n . 

Thus, a state on the composite system is denoted by a density operator p( n ) on HH In 
particular if n element systems (Hi, . . . ,H n ) of the composite system H^ n > are independent 
of each other, there exists a density p k on Hk such that 

pM = Pl ® •••<g>p n , onH {n) . 

The condition: 

Hi = ■ ■ ■ = H n = H, pi = ■ ■ ■ = p n = p 

corresponds to the independent and identically distributed condition in the classical case. 
Therefore, we consider the parameter estimation problem for the family {pj^ on H\6 £ 0} 
which is called the n-i.i.d. extended family. 

In this case, we use a sequence {M n } ( ^ ) =1 of POVMs where M n is a POVM on H whose 
measurable set is R d as an estimator. A sequence {M n }^ =1 is called an MSE consistent 
estimator if {M"}™ =1 satisfies 



lim^J \\0 -0\\ 2 trp e n) M n (d§) = 0, V# G 6. 



A recursive MSE consistent estimator can be regarded as an MSE consistent estimator be- 
cause a recursive estimator £ n = ({M k }^ =l , 6 n ) is regarded as a POVM M(£ n ) as follows: 

„ n 

M(£ n )(B):= I (g>M k (ui,...,u k „i)(dw k ), V£? C R d on H (n) . 

We define the quantum- correlational Cramer-Rao type bound Cg (G) for a weighted matrix 
G as: 

C® C (G) := inf {lminf nTr G%(M n ) {M"}™ =1 is an MSE consistent estimator } 

where the MSE matrix Vg(M n ) is given by: 

F> J '(M n ) = / (0* - - j ) tr pl n) M n ( d0). 

We have the following equation 

C® C (G) = lim inf nC$(G), 

where Cn(G) denotes the non-quantum-correlational Cramer-Rao type bound for the n- 
i.i.d. extended family [||. From the definition of the n-i.i.d. extended family, we have 
q Q y nCg(G). Therefore, we have the first inequality of ([!]). 

Cg QC (G) > C^ C {G) > C?(G). (1) 
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It shows the second inequality of ([[]) that Jg 1 " < nJg for any POVM M n on 7i^ n \ Therefore 
the difference between Cg QC (G) and C® C {G) means the difference of the quantum corre- 
lation case from the non-quantum correlation case. Under some regular condition, we can 
show that there exists an MSE consistent estimator {M n }^ =1 such that ||: 

nTiGV e (M n ) -> C® C {G) as n -» oo, V# 6 6. 



3 Quantum displaced thermal states family 

Now we consider the estimation for the the complex amplitude £ and expected photon number 
N for the quantum displaced thermal states family defined as: 



s := | PC ' N := / c exp (~^atM l«)(«l d2 « 



2' 



3.1 Estimation of complex amplitude £ 



C E C,JV > 



In the case of that photon number N is known, we estimate the tow unknown parameters 
( = (9 1 + i8 2 )/y/2. This estimation problem is investigated by Yuen, Lax and Holevo |^, 
In this case they calculated the inverse JT 1 of the RLD Fisher information matrix as: 



Jo 1 



They calculated the non-quantum-correlational Cramer- Rao type bound Cg (G) as follows: 



C^ QC (G) = Cf{G) = 2[N + ^g l + sjg\ - g\ - gl (2) 
where the weighted matrix G is parameterized as: 

q —( 51+52 53 \ 
\ 53 51 - 52 J ' 

From (||) and we have the following equations. 

Cq QC {G) = C® C {G) = C^{G). (3) 

In this case, the optimal estimator is the squeezed heterodyne. 

3.2 Simultaneous estimation of complex amplitude Q and ex- 
pected photon number N 

Next we consider the case of that both of the expected photon number N and the complex 
amplitude £ are unknown. In this case, we estimate three unknown parameters f? 1 ,^ 2 and 
# 3 = N. The first equation of (|3|) isn't held. Therefore, the squeezed heterodyne isn't 
optimal. The inverse Jg 1 of the RLD Fisher information matrix is calculates as: 

N + l § 
J^=\ -| N + ± 

N(N + 1) 
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Therefore we can calculated C^(G) as: 

C?(G) = g N(N + 1) + 2 (n + 51 + 1 <?! - 9 3 2 , 
if the weighted matrix G can be parameterized as: 

/ 91 + 92 93 \ 
G = 53 91 ~ 92 . (4) 
V 50 / 

If the weighted matrix G can be parameterized as (Q), we obtain the following equations: 

C% QC (G) >C$ C (G)=CK(G). (5) 

A proof for C^ QC {G) > C® C (G) is omitted. The inequality Cf QC (G) > C® C {G) means 
that we cannot the simultaneous measurement of the photon number counting and hetero- 
dyne for a single sample. 

3.3 Construction of an MSE consistent estimator {M n }™ =1 at- 
taining Cq-(I) 

Now, for a weighted matrix /, we construct an MSE consistent estimator {M n }^ =1 such that 

lim ntrV e (M n ) = C?(I), V0 G 9. 

n — >oo 

It is sufficient for C® C (I) = C^(J) to construct such an MSE consistent estimator. 
Every POVM M n is constructed in the following step: 

(1) Evolve the unknown state ® • • • (g> a s: 

V -— ■ v '-S 

n 

P(,N ® ■ ■ ■ ® PCAT PC,^ ® " " " ® PC,^ = Pv^C AT ® P0,Ar ® • • • ® /00,JV On ft (n) , 

v , ' " , ' " . ' 

n n n—l 

where 

U n = exp n _i(a* ai — a^a n ) • • • exp ^(a^oi — 0^03) exp ^(a^ai — 0*02) on W^ ra ^ 
1 

<pi = arctan — , i = 1, 2, . . . , n — 1. 

cii denotes the annihilation operator on Wj. 

(2) Measure the first sample Py^jy by the heterodyne, then we get the estimate £ of the 
complex amplitude. 

(3) Measure the others by the photon counting, then we obtain n—l data which obey 
the probability distribution P N (k): 

(4) We obtain the estimate ./V of the expected photon number N by the maximum likeli- 
hood estimator of the probability distribution P N (ki), . . . , P (k n -i). 
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4 Conclusion 



We formulate an asymptotic quantum estimation theory. This theory is applied to the 
simultaneous measurement of the photon number counting and the heterodyne for displaced 
thermal states. It is a future study to realize the MSE consistent estimator proposed in this 
paper in an actual physical system. 
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